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Abstract 

We study the cooling (heating) of a glue-parton gas due to production (destruction) of 
particles and determine the associated production of entropy. We incorporate sharing 
of the system energy among a changing number of particles. We find that the entropy 
of an evolving glue-parton gas changes in an insignificant range once the initial high 
temperature state has been formed, despite a great change in particle number and 
temperature. 
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In nuclear collisions carried out at energies per nucleon considerably greater than the 
nucleon rest mass we observe a very high final state particle multiplicity, with which is 
associated an unusually large entropy content Q . There are two questions which immediately 
come to mind when looking at this result: i) when and how is entropy produced in a quantum 
process, such as is a nuclear collision, and ii) how is the particle production related to entropy 
production? 

We will address here the second question and investigate in a quantitative manner the 
increase in entropy during the evolution of the primordial gluon-parton system 0, |3], |4| . In 
our approach, as in every prior work up to date, we do not view the nuclear collision reaction 
as a reversible quantum process, but rather as a series of disorder increasing independent 
reactions. Thus the measurement of (conventionally defined) entropy can be performed at 
any instant during the collisions |9, which approach may be seen in contradiction to the 
prevailing understanding of reversibility of quantum evolution. This is a pragmatic approach 
in view of the final state formed in the collisions, which is displaying high particle multiplicity. 
The view we and other workers take is that a solution to the more fundamental question i) 
will justify the pragmatic approach taken today regarding question ii), where we proceed as 
if the nuclear collision reaction was a classical process involving particle production. 

The physical systems of interest to us here are characterized by relatively large elastic and 
quasi elastic cross sections, which are typically 10-30 times larger than the particle number 
changing cross sections. Therefore it can be assumed that the so called kinetic equilibrium 
characterized by the sharing of the available energy between the particles, and the estab- 
lishment of a common temperature, is reached instantaneously as compared to the slowly 
evolving particle number. We call a particle species to be in Hagedorn equilibrium ||, when 
all phase space cells are occupied according to a distribution which maximize the entropy 
content at some given energy. It is interesting to note here that there is so far apparently 
no general proof of Boltzmann's H-theorem to be found in literature for relativistic systems 
with particle number non-conserving reactions ||. As expected we find below that Hage- 
dorn equilibrium corresponds to the standard Bose/Fermi quantum distribution of particles, 
which give a maximum of system entropy at fixed energy. This is a preliminary step towards 
such a proof, which requires further that the particle numbers, to be considered in kinetic 
theory, evolve in time towards the Hagedorn equilibrium — while for a select system |7|] this 
has been argued, a general discussion of this issue remains an open problem. In our detailed 
discussion we focus on the situation found in studies of the 'hot glue' initial state proposed 
by Shuryak 0. We will show that despite a great increase in particle number as the system 
approaches Hagedorn equilibrium relatively little additional entropy is generated. 

We begin with the considerations presented by Landau and Lifshitz ||. The usual defi- 
nition of entropy for each particle species denoted by index I is: 



where the integral is over the conventional phase space du = d^pd 3 x / (2nh) 3 (we henceforth 
choose units such that h — c — 1). The upper sign (here +) applies to the Bose particles 
('B'), while the lower sign (here — ) applies to the Fermi particles ('F'). The Boltzmann limit 
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follows when the occupation ni of the phase space by one of the particle species is small 
compared to unity. Given a number of particles of each species, N l = Jduni, and the energy 



of the system E = J2i E l = fdueini, where ei = ymf + p 2 , the form of ni(p,x) is obtained 
by maximizing the quantity T l [ni] = S l — (on + f3e{)ni as a functional of n\. a 1 permits 
to choose any given number of particles of kind I and (3 = 1/T, which permits to keep the 
energy constant, is identified as the inverse temperature. Reinserting the Bose 'B' and Fermi 
'F' quantum distribution functions 

which follow from this consideration into the definition of entropy, Eq. (|I|) we obtain the 
explicit form M: 



S 1 b,f —Jdu 



(3) 



eP e i+ a i =f 1 

The effect of a changing particle number appears in two different ways: as a factor 

li = e~ a < (4) 

in front of all Boltzmann exponential present in above equations and secondly, as a coefficient 
of an additional additive, particle number proportional term. We can thus obtain Si, Eq. (|3]), 
as a derivative of the partition function In Z ((3, o>i) = =f/ duj ln(l =p e~^ ei ~ a! ) . 

In the Boltzmann approximation the factor 7; becomes a normalization factor which 



describes the average occupancy of the phase space relative to the equilibrium value |TC| , the 
additive a^-term in Eq. (|3|) contributes to the entropy per particle change as the occupancy 
changes; in particular we note that in the Boltzmann limit: 

S l S l 



+wr 1 - ( 5 ) 



cq 



When we omit the last term from the discussion of the entropy of a system, we find that 
the entropy evolves along with the number of particles. Considering this term we see that 
as the occupancy of the particle phase space increases at fixed (3, the entropy per particle 
decreases. This observation is of course not in contradiction to our expectation that the 
total entropy should increase as the particle number increases. The situation appears to be, 
however, much less trivial when we look at what happens as the occupancy of phase space 
approaches the equilibrium value from above 1 < ji — > 1 . For the entropy must increase 
while the number of particles decreases. This can clearly only be the case if we consider 
how particle annihilation heats the system at given energy E. Similarly, we must consider 
when studying the increase of entropy during the approach from below to the equilibrium 
occupancy, how the production of particles cools the system. 

We find the required change in (3 for systems without conserved quantum numbers from: 

OE^ 1 OE 

dE = ^da- l dai+ dp dP ' (6) 
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which is inserted into the change in entropy to yield: 



JO JT ^dS/d0 _ J 

dS = dE MJdp + ^ dai 



dS l dSdE l /da t 
dai ~ df3 dE/dp 



(7) 



It is our hypothesis that the system is isolated, hence dE = 0. Considerable cancellation 
occurs when we evaluate Eq. (H): 



dS l 



duo 



-3d -/3e,EkJduj k a k e k -f k e 0ek ' 



one ^ ' — e;e 



J2 k JduJ k e k e k ^ k e-^ k 



(8) 



We see that indeed the change in entropy vanishes exactly when all oi\ vanish, that is 
when 7j — > 1. Furthermore we can see using Schwartz' inequality that for any 7; < 1 the 
change dS is positive for positive d'ji while for 7; > 1 the change is negative, proving that 
the full Hagedorn equilibrium (7; — > 1) gives maximum of entropy at fixed total energy. 

We consider now a systems of gluons in which numerous different reactions such as 
G* — ► GG — ► GGG^ • • • , are thought to increase the number of gluons rapidly and to 
lead to development of Hagedorn equilibrium. Our objective is to see the usually ignored 
effect of cooling due to particle production, and to determine the magnitude of the associated 
entropy increase. Our calculations are set up such that the Hagedorn equilibrium results at 
T = 0.250 GeV — due to energy conservation at fixed volume the energy density is constant 
as function of 7g (index 'G' is henceforth dropped) and takes the value 2.66 GeV/fm 3 . 
When we consider a thermal mass for gluons, which we take for illustrative purposes to be 
?71q = 0.200 GeV, we also choose to have Hagedorn equilibrium at T = 0.250 GeV and 
therefore the constant energy density drops to 1.89 GeV/fm 3 . 

We show in Fig. p] the entropy density S/V (units 1/fm 3 ) as function of 7. We note 
that for massless gluons we can not proceed beyond 7 = 1 because of Bose-condensation 
phenomena not considered further here. We show by dashed lines the same calculation 
carried out at = 0.200 GeV, which allows us to continue the calculation up to 7 ~ 2.7. 
The remarkable feature to be noted in the Fig. [I] is the appearance of the very weak maximum 
at 7 = 1 — we observe that for all practical purposes in the hot glue picture all entropy 
of the system is already created when the particle number evolution begins at 7 ~ 0.1 [|J. 
However, the temperature T shown in Fig.|2| evolves rapidly in the interval 0.1 < 7 < 1. 
While the total entropy increases by 15% (respectively 11%) the temperature drops from 
0.45 GeV (respectively 0.43 GeV) to 0.250 GeV for = 0. (respectively = 0.2 GeV), 
this decrease being driven here solely by particle production (no expansion cooling). 

As we can further note in Fig. ^| the glue number density increases considerably in the 
same region — it more than doubles for massless gluons. In Fig. |3] we also note that the glue 
number density increases monotonically as we pass 7 = 1 for itlq > 0. On the other hand, 
Fig. £| shows that the entropy per particle drops continuously as could be expected from the 
qualitative result obtained in the Boltzmann limit, Eq. (|5|). 

We note that in the expanding glue fireball there are no other known sources of entropy, in 
particular since the expansion is believed to be a largely entropy conserving hydrodynamical 
flow process. We will now address in qualitative terms the interplay between expansion and 
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Figure 1: Entropy density S/V (units 1/fm 3 ) 
at fixed energy density E/V = 2.66 GeV/fm 3 
for itiq = (solid line) and at E/V = 1.89 
GeV/fm 3 for m G = 0.200 GeV (dashed line) for 
a (gluon) Bose g clS ELS function of the chemical 
occupancy 7. 



Figure 2: The temperature T as function of 
the chemical occupancy 7. Lines as in Fig. |l|; 
equilibrium point 7=1 has been chosen to 
occur at T = 0.250 GeV. 



particle production cooling. In principle we need to know with some precision the relative 
rates of change of the volume V and the occupancy factor 7. Details of the such calculations 
are by necessity model dependent, but we can obtain some interesting qualitative insights 
pursuing some general relations. For (nearly) massless gas we have: 



S = - Cl 1/T 3 7 + c 2 1/T 3 7 a . 
3 



th 



En — E< 



flow 



Cl VT 4 ~f . 



(9) 



Both ci(7,/3m), and 02(7, (3m) are well known factors and can be exactly computed. For 
(3m — > we find: 



00 ^n-l 

y - — 

c 2 



00 n-1 




(10) 



n=l 



n" 



In the Boltzmann limit the factors C\, c-i are independent of a = In 7 1 . From Eq. @ we 
obtain using C\/oi ~ 3: 



S~S 



T 1 + q/ 4 
T 1 + oqj 



E{i ow \ 



/4 V E ) ■ 
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Figure 3: Particle density N/V (units 1/fm 3 ) 
as function of the chemical occupancy 7. Lines 
as in Fig. |]. 
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Figure 4: Entropy per particle S/N for a 
(gluon) Bose g clS cLS function of the chemical 
occupancy 7. Lines as in Fig. [[]. 



When -Eflow = this is a condition satisfied by the results we presented before. Because of 
the appearance of the flow factor in Eq. (jTT| ) which reduces the r.h.s. we further see that 
a reduction of the reservoir of thermal energy entails a further reduction in the growth in 
entropy due to particle production. 

Since entropy increase was rather small anyway, it is now possible to suppose that we 
have S/Sq — 1. Then we find from Eq. ( |Tl~D that the temperature cooling on the way to 
Hagedorn equilibrium is given by: 

T 1 + « /4 V E / 

where Th is the temperature at which the Hagedorn equilibrium is established, and there 
the fraction E^/Eq of energy was transferred to the flow. For an initial 70 = 0.1 occupancy, 
we have a = 2.3 and hence the cooling due to particle production is by factor 0.63 (correct 
numerical Bose gas result which accounts for entropy production is 0.55). The time period 
is of the order of 1 fm/c and during this time the fraction of energy transferred to flow may 
be as large as 40%. Consequently, we see that the particle production cooling (first factor 
in Eq. fll2|)) is of the same importance as the hydrodynamical expansion cooling (second 
factor in Eq. flT2|)) and that both combine to increase the temperature reduction, while both 
do not lead to appreciable entropy increase. We see that the cooling shown in Fig. is 
certainly of equal if not greater importance than the cooling expected from the longitudinal 
hydrodynamic expansion of the glue-parton volume obtained in a recent numerical study [|J . 

In summary, reviewing an old topic we have found that contrary to intuitive expectations 
the entropy growth is not proportional to the particle number increase, and it is indeed small 
during approach to Hagedorn equilibrium (7 — > 1) by a hot glue gas. We have used this 
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result to show that the temperature cooling due to particle production is comparable to 
cooling due to hydrodynamic expansion. We have seen in detail why a sizable increase in 
gluon number is not generating a large increase in entropy in the nuclear collision. Our 
finding thus renders the mechanism of entropy production in high energy nuclear collisions 
even less understandable, as it now must occur in the initial (quantum) time period before 
0.3 fm/c ||, when the primordial hot glue phase is generated. 
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